DESY 98-139 
|hep-ph/98094IB 

September 1998 



Double Spin Asymmetries in P-wave Charmonium 

Hadroproduction 



W.-D. NOWAK, A. TKABLADZEQ 

DESY Zeuthen, D-15738 Zeuthen, Germany 



Abstract 

We discuss the double spin asymmetries in P-wave charmonium hadroproduction with non- 
zero transverse momenta at fixed target energies, y/s ~ 40 GeV, within the framework of 
the factorization approach. The size of the asymmetries and the projected statistical errors 
in a future option of HERA with longitudinally polarized protons scattering off a polarized 
target (HERA-iV) are calculated. Measurements of the Xd and \c2 decays into dilepton 
plus photon should allow to distinguish between different parametrizations for the polarized 
gluon distribution in the proton. At higher energies (sfs = 200 GeV) the situation appears 
less favourable with the presently envisaged integrated luminosities of the polarized RHIC 
collider. 
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1 Introduction 

The study of spin asymmetries in the production of heavy quarkonium states in polarized 
nucleon-nucleon collisions should provide important information about the spin structure of 
the nucleon. Heavy quark-antiquark production processes occur at small distances and the 
subprocess level cross sections can be calculated perturbatively. On one hand, charmonium pro- 
duction asymmetries are expected to be sensitive to the polarized gluon distribution function in 
the proton, since heavy quark systems are mainly produced in gluon-gluon fusion subprocesses. 
On the other hand, it is essential to investigate in more detail the heavy quark-antiquark pair 
hadronization phase. To this end, observation of charmonium production in experiments with 
polarized beams is expected to provide additional tests for existing models. 

The Factorization Approach (FA) based on the Nonrelativistic QCD (NRQCD) turns out 
to be a most rigorous framework to study heavy quarkonium production and deacy processes 
0. According to the FA the inclusive production cross section for a quarkonium state H in the 
process 

A + B -> H + X (1) 

can be factorized as 

a(A + B^H)=j:^(0\O»\0), (2) 

n m Q 

where the short- distance coefficients, F n , are associated with the production of the heavy quark 
pair in the color and angular momentum state [n]. This part of the production cross section 
involves only momenta of the order of the heavy quark mass and larger and can be calculated 
perturbatively. The heavy quark-antiquark pair production process occurs at small distances, 
(1/mg), and is factorized from the hadronization phase which takes place at large distances, 
1 / (itiqv 2 ), where uiq is the heavy quark mass. Here v is the average velocity of heavy constituents 
in the quarkonium, with v 2 ~ 0.3 for charmonium and v 2 ~ 0.1 for bottomonium systems. The 
vacuum matrix elements of NRQCD operators, (0|C^|0) JTJ, describe the evolution of the quark- 
antiquark state [n] into the final hadronic state H. These matrix elements cannot be calculated 
perturbatively, but the relative importance of long distance matrix elements in powers of velocity 
v can be estimated using the NRQCD velocity scaling rules [J2J. This formalism implies that 
quark-antiquark color octet intermediate states are allowed to contribute to heavy quarkonium 
production and decay processes at higher order in the velocity expansion. Therefore, in the 
FA the complete structure of the quarkonium Fock space is taken into account while in the old 
approach, in the Color Singlet Model (CSM) ||, only the dominant Fock state is considered, 
which consists of a color singlet quark-antiquark pair in a definite angular-momentum state of 
the final hadron (the leading term in the velocity expansion). 

The predicted shape of the pr distribution of the 3 Sf^ intermediate octet state production 
cross section indicates that J/ ip and ■0' production at large p T observed at the Tevatron (FNAL) 
can be explained in the FA || ||. Recent investigations have shown that the contribution of 
color octet states to the charmonium production cross section is very important at fixed target 
energies, yfs ~ 30 — 60 GeV, and reduces existing discrepancies between experimental data and 
predictions of the CSM g. 

Despite the obvious successes of the NRQCD in explaining large pt charmonium production 
at the Tevatron, some experimental data contradict the Color Octet Model (COM) predictions. 
In particular, theoretical predictions disagree with measurements of the polarization of J/ip and 
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if)' particles produced at fixed target energies || and the COM prediction for the yield ratio 
of Xci an d Xc2 states remains too low ||. One possible solution for these discrepancies was 
proposed by Brodsky et al. |J suggesting that higher twist processes, when more than one 
parton from projectile or target participate in the reaction, might give a significant contribution 
to low pt production of J / if) and Xci states. Problems exist also in charmonium photoproduction 
at HERA. The color octet contribution underestimates the inelastic J /if) photoproduction cross 
section at large values of z (z — Ej/^/E^ in the laboratory frame) [|TC|j . 

Unlike the color singlet long distance matrix elements, each connected with the subsequent 
hadronic non-relativistic wave function at the origin, color octet long distance matrix elements 
are unknown and have to be extracted from experimental data. The NRQCD factorization 
approach implies universality, i.e. the values of long distance matrix elements extracted from 
different experimental data sets must be the same. However, due to the presently rather large 
theoretical uncertainties JO], |T3| and the unknown size of higher twist process contributions || 



the existing experimental data does not allow to check the FA universality, yet. 

This fact motivated us to look for other processes with less theoretical uncertainties to test 
the COM. The observation of spin asymmetries in the production of charmonium states can 
be used for these purposes |L2| as well as measurements of the J /if) polarization in unpolarized 



hadron-hadron collisions and electroproduction 0, O, P 



In the present article we consider double spin asymmetries in P-wave charmonium production 
at non-zero pt in the NRQCD factorization approach. The double spin asymmetry in J /if) 
hadroproduction has been studied recently taking into account color octet intermediate states 
|T2| , [R| . The asymmetries in the production of Xcj states at small px were considered in refs. 



15|, [Tq| . As already mentioned above, at small values of px large contributions from higher twist 
effects are expected and the presently available theoretical predictions are not reliable enough 
to allow extraction of information about polarized parton distributions or to check the COM. 
Morii et al. considered asymmetries in the production of Xcj states at RHIC energies taking into 
account only fragmentation type contributions for color octet intermediate states (I7j . However, 



since the cross section for the production of charmonium states falls steeply with increasing p? 
the statistical errors becomes too large at those values of pt where the fragmentation approach 
can be applied, especially at at small energies ^fs ~ 50 GeV. In contrast to this approach, our 
calculations |H| use the exact expressions for the cross sections of intermediate color octet states 
in different helicity states of the initial partons, to estimate the expected spin asymmetries in 
the production of Xcj states as well as the projected errors. 

The article is organized as follows. In the next section subprocess level asymmetries for possi- 
ble intermediate octet and singlet states and corresponding long-distance color octet parameters 
are considered. In section 3 numerical results are presented for spin asymmetries in production 



of XcJ states at HERA-iV ||23|| , one of the future options of HERA; an experiment using an 
internal polarized nucleon target in a possibly later polarized HERA beam at ^Js ~ 40 GeV. 
For comparison, we also consider the expected asymmetries in the production of Xcj states in 
similar spin physics experiments at the RHIC collider at -^i = 200 GeV. 

2 Asymmetries at the Subprocess Level 

The two-spin asymmetry All for inclusive XcJ production is defined as 

A LL (pp) = da(j)+P+ — XcJ ^ - da ( p+p - Zl XcJ "> = EdAa / d3 P / 3 x 
LL dcr(p + p + ^Xcj) + da{p+p- -> Xcj) Eda/d 3 p 



3 



where the subscript in indicates the sign of the helicity projection onto the direction of 

the proton. The production of each quarkonium state receives contributions from both color 
octet and color singlet states. In leading order v 2 only one color octet state contributes to the 
production of Xcj states, namely 3 S[ 8 \ Unlike direct J/ip production, leading color octet and 
color singlet contributions to P-wave charmonium production scale equally in v 2 , 0(v 5 ), and the 
subleading corrections are only of the order 0(v 9 ). To calculate the production asymmetries 
only leading order color octet contribution are taken into account. 





Figure 1: The subprocess level asymmetries for 3 Pq, 3 Pi, and 3 P2 color singlet, as well as for 3 S[ 
color octet state. 

We consider the production of x states at non-zero values of transverse momentum, > 1.5 
GeV. The leading contribution to the production of charmonium states at those values of trans- 
verse momentum comes from the 2^2 subprocesses. The leading order subprocesses, 2 — > 1, 
should not contribute to the quarkonium production with p T > 1.5 GeV because these transverse 
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momenta cannot be caused by internal motion of partons in the nucleon. The contribution of 
higher twist processes M is also expected to be small at px > 1.5 GeV. 

For the calculation of the expected asymmetries for color octet and singlet states of a (cc)-pair 
we consider the following subprocesses: 



9 + 9 

g + q 

q + q 



{cc) + g 
{cc) + q 
{cc) + g 



(4) 



The asymmetries for all these subprocesses were calculated in ref . |12| . The values of those giving 
a significant contribution to the production of Xcj states are presented in Fig. 1 in dependence 
on the two dimensionless quantities rj = Am 2 / s and xt = Pr/Pmax, where m c is the charm quark 
mass, s the cm. energy in the parton-parton system, and p ma x is the maximum momentum of 
the produced state in the subprocess. Figure 1 shows only the gluon-gluon fusion subprocess 
asymmetries because they give the main contribution to the hadron level asymmetry. 

The color singlet long distance matrix elements are connected to quarkonia wave functions 
at the origin: 



(o\or j ( 3 Pj)\o) 



2vr 



(2J + l)|i?'(0)| 2 = 3.2 ■ lO^GeV 3 . 



(5) 



This value was used in J 



for fitting the CDF data on J ftp production through Xcj states and 

For the color octet 



corresponds to the Buchmuller-Tye type potential solution tabulated in |21 



long distance matrix element the following value extracted from CDF data was used 

(OlOf^SOlO) = 9.8 • 1(T 3 GW 3 , 



(6) 



0.8 



xAG(x,Q 2 =10GeV 2 ) 




Figure 2: Different possible polarized gluon 
distributions in the nucleon, used in this pa- 
per: GS parametrization NLO set A (solid 
line), NLO set B (dashed line) and LO set A 
(dash-dotted line) §§. 



As shown in ref. [O] the variation of the renormal- 
ization and/or the factorization should lead to large 
uncertainties when fitting the color octet parameters. 
The fit results for long distance color octet matrix 
elements can be affected also by higher v 2 correc- 
tions (so called 'trigger bias' effect) [|TTJ. Therefore 
care is required when using at fixed target energies 
the value for this parameter extracted from CDF. 
The influence of these uncertainties on the expected 
asymmetries will be discussed in the next section. 

3 Results and Discussion 

The characteristic value of the partonic x in the pro- 
duction of (cc) pairs can be obtained from the re- 
lation x x x 2 ~ (4m 2 +p 2 T )/s (~ 0.01 at HERA-AO. 
This means that the typical value of x g i uon that can 
be probed by measuring the spin asymmetry in char- 
monium production is about Xgi uon ~ 0.1. We used 
three different sets of Gehrmann and Stirling (GS) 
parametrizations |22| for polarized parton distribu- 
tion functions (PDF) that are different in the region 
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x ~ 0.1 to show the sensitivity of the spin asymmetries in the production of Xc.j states on the 
gluon polarization in the nucleon: the sets A and B of the NLO GS parameterization and the 
LO parameterization set A. In Fig. 2 the polarized gluon densities from these parametrizations 
are shown at Q 2 = 10 GeV 2 , which is the appropriate value in the kinematical region considered. 
As can be seen, the three chosen sets exhibit somewhat different values for the polarized gluon 
distribution function at partonic x values around 0.1. We note that, although the calculations of 
subprocess level cross sections are performed in leading order, the NLO set of the parametrization 
was used to probe different shapes of the polarized gluon distribution. 




p T (GeV) p T (GeV) 



Figure 3: a) Xd an d b) Xc2 production asymmetries versus transverse momentum at y/s = 40 GeV. The 
solid (dashed) line correspond to set A (B) of the NLO GS parametrization and the dash-dotted line 
to set A of the leading order GS parametrization. The dotted line represents the expected asymmetry 
calculated only in the CSM for the set A of the NLO parametrization. 

Figs. 3 and 4 show the expected production asymmetries for Xcj states at y/s = 40 GeV 
(HERA-iV). Figures 3a and 3b also include projected statistical uncertainties on the asymmetry 
5A LL , estimated at HERA-iV from @ by 

6A LL = 0.17/ sj^pb)- (7) 

where 100% efficiency is assumed. This relation is based upon an integrated luminosity of 
240 and beam and target polarizations of Pb = 0.6 and Pt = 0.8, respectively p3| . The 
error bars are obtained by using integrated cross sections over bins Apx = 0.5 GeV (for the 
first three points) and Ap? = 1 GeV (for the other two ones). Note that it is easier to measure 
asymmetries for Xci an d Xc2 states due to their large branching ratios for decays into J /if) plus 
photon (27% for Xic and 13.5% for X2c)- These branching ratios and the J /if) decay branching 
ratios into (I = e, fi) are taken into account in the calculations of the projected statistical 
errors for the production of Xci and Xc2 states. Due to the small branching ratio of the Xco decay 
through J/if>, 6.6 • 10~ 3 , the projected statistical errors of an asymmetry measurement for this 
state are too large. 

In Figs. 3 and 4, the asymmetries for set A and set B of the NLO GS parametrization are 
represented by solid and dashed lines, respectively, whereas the dash-dotted lines correspond to 
set A of the LO parametrization. The asymmetries coming only from color singlet intermediate 
states expected for the set A of NLO GS parametrization are additionally shown in Figs. 3 and 
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4 by dotted lines. Note that for the set C of the NLO GS parametrization all asymmetries are 
less than 1% and practically unobservable, hence they are not shown. For the mass of the charm 
quark the value m c = 1.48 GeV was taken and the parton distribution functions are evaluated 
on the factorization scale \x - 
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lp\ + Am? c . The strong coupling constant is calculated by the 
one-loop formula with 4 active flavors (Aqcd = 200 MeV). 

In all three parameterizations for the polarized parton distribution functions the gluon-gluon 
fusion process gives the dominant contribution to Act in the kinematical region considered. The 
quark-gluon subprocesses contribute only to about 10% (for both color octet and color singlet 
states) and the contribution of quark-antiquark annihilation subprocesses is less than 1%. 

The contribution of the color octet 3 Sj 8 ' ) state 
does practically not change the Xci production asym- 
metry coming from the color singlet intermediate 
state. For Xco, the influence of the octet contribu- 
tion becomes larger with increasing px- As seen from 
o.i5 h H Fig. 1, the asymmetries for color singlet 3 P2 and color 

octet 3 S[ 8 ^ states have different sign. By this reason 
the color octet contribution, compared to the color 
o.i h H singlet one, changes the asymmetry of the Xc2 state 

more than in the case of Xco and Xci- Moreover, due 
to the 'trigger bias' effect the corresponding color 
o.os h H octet long distance matrix element may be approxi- 

mately two times larger and hence the contribution 
of the color octet states to the asymmetries may even 
increase. 

In the inclusive case the kinematics of the 
2 — ^ 2 subprocess cannot be reconstructed com- 
pletely. Hence only indirect information on the gluon 
polarization of the nucleon can be obtained by mea- 
suring the spin asymmetry in inclusive Xcj produc- 
tion. As can be seen from Figs. 3 and 4, the expected 
asymmetries for all Xc states strongly depend on the size of the polarized gluon distribution func- 
tion in the region x g \ uon > 0.1. At HERA-iV the size of the Xci and Xc2 asymmetries compared to 
the projected statistical errors will allow to draw conclusions on the polarized gluon distribution 
function. 

In addition, the measurement of xa and Xc2 asymmetries provides a good opportunity to 
discriminate between two possible mechanisms of heavy quarkonium production: the NRQCD 
factorization approach and the Color Evaporation Model (CEM) ||24j| . As can be seen from 
Fig. 3b, the Xc2 production asymmetry is negative at smaller values of transverse momentum 
where the gluon-gluon fusion subprocess contribution to Act is dominant. In spite of the positive 
contribution of the color octet state 3 S[ 8 ^ to Act, the Xc2 production asymmetry remains negative 
as in the CSM, i.e. in the NRQCD approach the double spin asymmetries for Xci and Xc2 states 
have different sign. In contrast, in the CEM the asymmetries for all charmonium states are 
expected to be the same. Such a comparison is possible even when shape and size of the 
polarized gluon distribution in the proton are not yet known. The only requirement is that the 
gluon polarization is large enough to generate observable xa and Xc2 asymmetries. 

For comparison we have also calculated the expected double-spin asymmetries for the produc- 
tion of Xci and Xc2 states at RHIC energies. The results are given in Fig. 5 for the c.m.s. energy 



4 4.5 5 
p T (GeV) 



Figure 4: Expected asymmetries for Xco pro- 
duction at HERA-iV for different GS param- 
eterizations as explained in Fig. 3. 
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*/s = 200 GeV. The statistical errors are calculated with the anticipated integrated luminosity 
of 320 pb~ x assumming 100% efficiency and PbPt — 0.5. As can be seen by comparing Fig. 3 
and Fig. 5, the expected asymmetries decrease with increasing c.m.s. energy for a given pa- 
rameterization of polarized parton distribution functions. Moreover, the anticipated integrated 
luminosities at HERA-iV and RHIC lead to a smaller discrimination power at higher energy. 




p T (GeV) p T (GeV) 

Figure 5: a) Xci and b) Xc2 production asymmetries versus transverse momentum at = 200 GeV. 
The solid (dashed) line correspond to set A (B) of the NLO GS parametrization. The dotted line 
represents the expected asymmetry calculated only in CSM for the set A of the NLO parametrization. 



One of the main parameters of the factorization approach is the mass of the charm quark. 
As in the case of J/ip production fl2| , the expected asymmetries for Xcj states are practically 
insensitive to the value of the charm quark mass. Therefore, the double spin asymmetry in Xcj 
production, unlike the cross section, should be free from uncertainties caused by the unknown 
mass of intermediate color octet states. Moreover, the asymmetries do not strongly depend 
on the renormalization scale, unlike the production cross sections [|T3J. The effects of intrinsic 
transverse momentum smearing on the asymmetries in the production of charmonium states were 
also investigated. Such an effect is very important for the cross section of J/ip photoproduction 



at HERA collider energies [25], however its influence to asymmetries is far below the size of the 



projected statistical errors shown in Figs. 3 and 5. 



4 Conclusions 

Double spin asymmetries expected in hadroproduction of heavy quarkonium P-wave states in 
polarized proton-proton collisions were investigated. To reduce the contribution from possible 
higher twist corrections M the production of Xcj mesons was considred at non-zero transverse 
momenta, pr > 1-5 GeV unlike the calculations of | 15| , |I6|, where only the lowest order subpro- 
cesses 2 — > 1 were taken into account. 

The size of the expected asymmetries in conjunction with the projected statistical uncertain- 
ties at HERA- A" (-y/s ~ 40 GeV) will allow to distinguish between different parametrizations for 
polarized parton distribution functions, specifically between the NLO GS sets A and B [p^j . At 
RHIC, higher than hitherto planned integrated luminosities would be required to do so. 

The asymmetries for the production of the xa and Xc2 states in the transverse momentum 
range 1.5 < pt < 3 GeV exhibit different signs in the NRQCD factorization approach. In 
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contrast to this, the color evaporation model predicts identical double spin asymmetries for all 
charmonium states. This should allow to discriminate between these different mechanisms of 
heavy quarkonia production if the gluon polarization is sufficiently large and Xc asymmetries 
become observable. 

We are grateful to Oleg Teryaev for helpful discussions. We thank J. Bliimlein for careful 
reading of the manuscript. A.T. acknowledges the support by the Alexander von Humboldt 
Foundation. 
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